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Abstract 

Purpose: The mechanism underlying tlie decline in milk quality during periods of feeding high-concentrate diets to dairy 
ruminants is not well documented. The aim of this study was to investigate the metabolic changes in the liver that 
contribute to the input of substrate precursors to the mammary gland after feeding a high-concentrate diet to lactating 
goats for a long period. 

Experimental Design: Eight mid-lactating goats with rumen cannulas were randomly assigned to two groups. For 9 weeks, 
the treatment group was fed a high-concentrate diet (60% concentrate of dry matter, HC) and the control group was fed a 
low-concentrate diet (40% concentrate of dry matter, LC). Ruminal fluid, plasma, and liver tissues were sampled, microarray 
techniques and real-time polymerase chain reaction were used to evaluate metabolic parameters and gene expression in 
liver. 

Results: Feeding a 60%-concentrate diet for 9 weeks resulted in a significant decrease in rumen pH. Changes in fat and 
protein content also occurred, which negatively affected milk quality. Plasma levels of leptin (p = 0.058), non-esterified fatty 
acid ip = 0.071), and glucose (p = 0.014) increased markedly in HC group. Plasma Cortisol concentration was significantly 
elevated in the treatment group {p<0.05). Expression of the glucocorticoid receptor protein gene was significantly down- 
regulated {p<0.05) in the liver. The expression of genes for interleukin 1 p, serum amyloid A, C-reactive protein, and 
haptoglobin mRNA was significantly increased (p<0.05) in the HC group. GeneRelNet analysis showed that gene expression 
involved in inflammatory responses and the metabolism of lipids, protein, and carbohydrate were significantly altered by 
feeding a high-concentrate diet for 9 weeks. 

Co/Jtr/i/5/b/75; Activation of the acute phase response and the inflammatory response may contribute to nutrient partitioning 
and re-distribution of energy in the liver, and ultimately lead to a decline in milk quality. 
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Introduction 

Dairy ruminants are often fed a liigh-concentrate diet to meet 
the energy requirements of high milk production. However, 
metabolic disorders such as subacute ruminal acidosis (SARA) can 
result from feeding excessive amounts of high proportions of 
fermentable concentrate and forage with low physically effective 
fiber [1]. SARA is a common digestive disorder in most dairy 
herds. In intensive management systems, 19% of early lactation 
and 26% of mid-lactation cows can be affected by SARA [2]. The 
consequences of SARA include depression of food intake, 
diarrhea, laminitis, and inflammatory response, which ultimately 
result in the depression of milk quality and quantity [3] [4] [5] [6] . 



Production costs from SARA have been estimated to be as high as 
$1.12 per cow per day [7] . SARA is a major concern for the global 
dairy industry. Especially, in the regions shortage of high quality 
forage and lack of reasonable feeding management practices 
prevents dairy cows from achieving their genetic potential for milk 
production. 

SARA is experimentally characterized by repeated bouts of 
depressed ruminal pH values between 5.2 and 5.6 for at least 
3 hours per day or below5.8 [5] [6]. Overfeeding fast fermentable 
starch or feeding excessive concentrate to forage ratios for 2 to 
3 weeks can induce SARA clinical signs [2]. Most studies have 
focused on the short-term (2-3 weeks) effects of SARA on 
inflammatory response, rumen fermentation, and milk production. 
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Table 1. Comparison of plasma concentrations of glucose, 
non-esterified fatty acids (NEFA), triglyceride (TG) and leptin. 



Diet 



Parameters 


LC 


HC 


/'-value 


Glucose 


3.00±0.06 


3.3 ±0.05 


0.014 


NEFA (nmol/L) 


175.87±55.12 


969.52±293.20 


0.071 


TG (mmol/L) 


0.19±0.06 


0.1 8± 0.03 


0.659 


Leptin (ng/mL) 


T75±0.57 


2.01 ±0.78 


0.058 



LC, low concentrate; HC, high concentrate; NEFA, non-esterified fatty acid. TG, 
triglyceride. Values are mean ± standard error of the mean (SEM), n = 4/group. 
doi:l 0.1 371 /journal.pone.0082850.ta01 



Fairfield et al. [7] and Kliafipoor [8] reported that experimentally 
induced SARA decreased milk fat percentage but increased milk 
protein content. However, other reports have suggested that 
during the short term SARA did not affect milk fat or protein 
content [9] [10]. However, dairy production herds affected by 
SARA can experience negative outcomes as long as several weeks 
(>3 weeks) after an episode [1 1]. These diseases include laminitis, 
weight loss, and unexplained abscesses. The long-term metabolic 
effects of SARA have not been well studied. 

Feeding high-concentrate diets to lactating cows results in the 
release of bacterial endotoxins, such as lipopolysaccharide (LPS), 
from the rumen or hind gut [12]. Free LPS can translocate into 
the bloodstream from the digestive tract under conditions of high 
permeability and after injury to the epithelial tissue [12]. Immune 



responses are then activated and pro-inflammatory cytokines are 
released from toll-like receptors 4 (TLR-4) that recognize the 
circulating LPS [13]. The process of LPS recognition via TLR-4 is 
facilitated by the accessory molecule LPS-binding protein (LBP) 
and cluster of differentiation antigen 14 (CD 14) [14] [15]. LBP is 
an acute phase protein (APP). The acute phase response to 
infection, tissue injuries, or other disorders participates in the 
restoration of physiological homeostasis [16]. The liver is the 
major site for the synthesis of APPs, including LBP, C-reactive 
protein (CRP), serum amyloid A (SAA), and haptoglobin (HP) [17] 
[18] [19]. Monokines such as interleukin 6 (IL-6), IL-1, and tumor 
necrosis factor alpha (TNF-a) induce APP synthesis via their 
specific hepatic receptors [20] These monokines also stimulate tiie 
release of adrenocorticotropic hormone (ACTH) from the 
pituitary gland. ACTH signals an increased secretion of gluco- 
corticoid hormone from the adrenal cortex [21]. During the 
restoration of physiological homeostasis, these mediators act on 
specific receptors on different target cells, which results in changes 
in metabolism and systemic responses. 

As the vital organ that controls metabolism [22], the liver is 
stimulated by endogenous and exogenous factors. It is responsible 
for nutrient partitioning, including homeostatic and long-term 
homeorhetic adaptations [23]. Physical and psychological stress 
elevates plasma IL-6 and APP levels in humans and in 
experimental animals [24] [25]. Although the mechanism of 
APP induction in response to stress is yet to be elucidated, 
activation of the hypothalamic-pituitary-adrenal (HPA) axis by 
stressors may trigger systemic or local cytokine production. It's 
reported that, during stress response physiological processes are 
aimed on redistribution of energy utilization in specific organs. 
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Figure 1. Rumen fluid pH values. Data were analyzed for differences due to diet, time, and their interactions by Univariate using the General 
Linear IVIodels of SPSS 11. 0 for Windows (StatSoft Inc, Tulsa, OK, USA). Values are mean ± standard error of the mean (SEM). 
doi:1 0.1 371/journal.pone.0082850.g001 
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Table 2. The percentage of fat, protein, and lactose in milk from control and treatment group goats. 



Item Diet ''"ime Effect, P-value 

0-2 W 3-4 W 5-6 W 7-9 W Diet Time Diet x Time 

Fat LC 4.15±0.43 4.01 ±0.38 3.76±0.88 3.90±0.53 

HC 3.09±0.63* 3.34±1.15 3.20±1.22» 3.04± 0.89** 0.931 0.021 0.935 

Protein LC 2.90±0.35 3.26±0.16 3.10±0.67 3.18±0.10 

HC 2.88±0.28 2.92±0.40 2.62±0.4n 2.68±0.22** 0.516 0.007 0.372 

Lactose LC 4.82±0.23 4.59±0.33 4.72±0.28 4.78±0.29 

HC 4.90±0.28 4.96±0.09« 5.00±0.10 4.95±0.13* 0.845 0.017 0.646 

LC, low concentrate; HC, high concentrate; Values are mean ± standard error of the mean (SEM), n = 4/group. * p<0.05, ** P<0.01 vs. LC. 
doi:1 0.1 371 /journal.pone.0082850.t002 



inhibiting or stimulating energy mobilization according to priority 
[26]. As a fact, defense against infectious diseases and stress 
response are resource expensive [27]. Tlierefore, the objective of 
this study was to investigate the changes of acute phase response, 
stress and inflammatory response, as well as the metabolic changes 
in liver that may contribute to the input of substrate precursors to 
the mammary gland after feeding a high-concentrate diet to 
lactating goats for a long period. 

Results 

Ruminal fluid pH 

After 1 week of feeding a high-concentrate diet, the mean 
ruminal pH value in the HC group was pH <5.8 from 1 to 
5 hours after feeding, which was consistent with the experimental 
definition of SARA (figure lA). During week 4 to week 8, enhanced 
ruminal buffering capacity and diet adaptation resulted in an 
increase in pH of HC group (mean pH >6.0), but it was still lower 
than that of LC group (figure IB, C and D). 

Plasma biochemical parameters and milk composition 

In HC goats compared with LC counterparts, the levels of 
plasma non-esterrfied fatty acid (NEFA; p = 0.071) and leptin 
(J) = 0.058) had a tendency to increase and blood glucose 



Table 3. Lipopolysaccharide (LPS) levels in plasma and in 
rumen fluid of control and treatment group goats. 







Items Diet 




SEM Effect, /'-value 










Diet X 


LC 


HC 


Diet Time 


Time 


Plasma LPS, EU/mL 


Oh 0.24 


0.32 


0.05 0.249 0.002 


0.716 


4h 0.09 


0.14 






8h 0.07 


0.06 






Rumen LPS', EU/mL 


Oh 38,854 


48,064 


20,383 0.717 0.084 


0.98 


4h 14,466 


16,058 






8h 63,359 


70,418 









LC, low concentrate; HC, high concentrate; LPS, lipopolysaccharide; Values are 
mean ± standard error of the mean between two groups (SEM), n = 4/group. 
'statistical analysis conducted on logiQ— transformd date. 
doi:l 0.1 371/journal.pone.0082850.t003 



{p = 0.014) was significantly increased after feeding high concen- 
trate diet for 9 weeks (Table 1). 

The percentage of milk fat and protein decreased significandy 
from 5-9 weeks in HC group compared with LC group {p<0.05). 
However, mUk lactose concentration increased significantly in HC 
group from the 4th week of treatment compared to LC group 
(p<0.05; Table 2). 

LPS content in ruminal fluid and plasma 

LPS levels in ruminal fluids of HC group increased from 38,854 
to 48,064 EU/mL before feeding in the morning and increased 
from 63,359 to 70,418 EU/mL after 8 h feeding compared to LC 
group. Plasma LPS concentration of HC goats was increased from 
0.24 to 0.32 EU/ mL before feeding in the morning and from 0.09 
to 0.14 EU/mL after 4h feeding compared with LC group. 
However, there was no significant difference in plasma LPS 
concentration and ruminal LPS content between HC and LC 
group (p>0.05; Table 3). 

Cytokines, acute phase proteins 

Plasma concentrations of IL-6 and TNF-ot increased in HC 
group but were not significantly different compared to the LC 
group {p>0.05; Table 4). Compared with LC group, plasma SAA 
concentration in HC group was increased significantly {p<0.05), 
and plasma HP concentration also showed a tendency to increase 
{p = 0.08). Results for plasma LBP were not significantly different 
for the HC group compared to the LC group {p = 0.163). 



Table 4. Plasma cytokines and acute phase protein levels in 
control and treatment group goats. 



Diet 


Parameters 


LC 


HC 


/'-value 


IL-6(|.ig/mL) 


221.32±42.63 


345.45±56.33 


0.163 


TNF-0!{ng/mL) 


715.42±49.67 


984.36 ±130.46 


0.126 


SAA(ng/mL) 


90.62± 15.67 


246.04±37.37 


0.020 


HP{^g/mL) 


209.50±21.99 


307.54±35.77 


0.087 


LBP(ng/mL) 


23.32±6.10 


34.48±5.74 


0.163 



LC, low concentrate; HC, high concentrate; HP, haptoglobin; IL-6, interleukin-6; 
LBP, lipopolysaccharide binding protein; SAA, serum amyloid A protein; TNF-ct, 
tumour necrosises factor a. Values are mean ± standard error of the mean 
(SEM), n = 4/group. 

doi:l 0.1 371 /journa].pone.0082850.t004 
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Figure 2. Plasma Cortisol concentration (A) and hepatic 
glucocorticoid receptor (GR) protein expression (B). The 

concentrations of Cortisol were significantly up-regulation and the GR 
protein expression was significantly down-regulation. Values are mean 
± standard error of the mean (SEM); *p<0.05 vs. LC. 
doi:1 0.1 371 /journal.pone.0082850.g002 

Cortisol and glucocorticoid receptor 

The level of plasma Cortisol was significantly elevated in HC 
group (p<0.05). Glucocorticoid receptor (GR) protein expression 
in the liver was greatly down-regulated in HC group compared 
with LC group {p< 0.05; figure 2). 

Liver mRNA Expression Profiling by Microarray 

Approximately 512 genes were observed to be diflFerentiaUy 
expressed in the HC group with approximately two-thirds of the 
genes exhibiting increased expression in liver. After removal of 
redundant and unannotated sequences {p<0.01), 43 genes were 
significantly down-regulated and 1 1 2 genes were significantly up- 
regulated in treated goats compared to LC group {figure 3). 
Functional annotation of enriched Gene Ontology terms indicates 
that a number of biological processes within the liver are affected 
by consumption of high-concentrate diets, including those related 
to control of cell surface receptor linked signal-transduction, G- 
protein coupled receptor protein signaling, regulation of tran- 
scription, immune response, and proteolysis. 

The raw data and processed files for the two libraries have been 
deposited in NCBLs Gene Expression Omnibus database (http:// 
www.ncbi.nlm.nih.gov/geo/) under accession number GPL17623 
and GPL17624. 



The GeneRelNet (Co-expression network) analysis 

Co-expression network analysis was used to fmd K-core 
regulatory factors (genes) (figure 4A-D). For each pair of genes in 
a cluster, we calculated the Pearson correlation coefficient and 
chose the significant correlation pairs to construct the network 
(absolute value of interaction >0.999). These genes were mainly 
involved in immune and inflammatory responses (n = 8), lipid 
metabolism (n = 5), protein metabolism (n = 26), carbohydrate 
metabolism (n = 7), out of a total of 155 differential expression 
genes {p<0.01). 

Real-time polymerase chain reaction (PCR) analysis 
involved in inflammatory response, glucose metabolism 
and lipid metabolism 

The mRNA expression of hepatic TLR-4, IL-ip, SAA, CRP, 
and HP were significantly increased in HC group compared to LC 
counterparts (p<0. 05; figure 5A). The expressions of hepatic TLR2, 
CD 14, MyD88, TNF-a, and LBP did not show significant 
difference between HC and LC group {p>0.05). Among these 
genes, only HP mRNA expression showed a more than 2-fold 
difference between the two groups, which was consistent with the 
data from the transcriptome microarray. 

Expression of hepatic phosphoenolpyruvate carboxykinase (PC) 
mRNA was significantly down-regulated in HC group compared 
to LC group (p<0.05) figure 5B). Hexokinase (HK) mRNA 
expression had a tendency to decrease in HC group (p = 0.07). 
However, other gluconeogenesis- and glucolysis-related genes 
transcription was not significantly altered in HC group compared 
to LC group {p>0.05). 

The mRNA expression of hepatic acetyl-CoA acyltransferase 2 
(ACAA2), acyl-CoA synthetase long-chain family member- 1 
(ACSLl), and hydroxyacyl-CoA dehydrogenase (HADH) genes 
were significantly increased (p<0.05) in HC group compared to 
LC group figure 5Cj. There was a >2-fold difference in hepatic 
fatty acid desaturase 1 (FADSl) mRNA expression, but did not 
reach statistical significance due to the high variability between 
individual values. This result was consistent with the data from the 
microarray high-throughput data analysis. 

Discussion 

Dynamic changes in ruminal pH during SARA 

Dairy ruminants are often fed high grain diets to meet the 
energy demand for high milk production. As a result, ruminal 
acidosis (especially SARA) often occurs during the periods of early 
and mid-lactation in dairy production herds [13] [11]. Due to the 
rapid fermentation and the accumulation of volatile fatty acid 
(VFA) in the rumen, as well as the increase of lactic acid, the value 
of rumen pH markedly decreased and prolonged for long time 
[28] [29]. In experimental study, rumen pH value has been 
considered as the most effective indicator of SARA. The 
thresholds of ruminal pH for the diagnosis of SARA have been 
set as ruminal pH <5.6 [10] [30] [31] and <5.8 [28] lasted for 
3 hours per day using ruminal fluids detected by pH meter system 
in dairy cows [10]. Consistently, ruminal pH value <5.8 lasted 
more than 3 hours was considered as the critical value of SARA 
diagnosis [32]. 

In this study, goats fed a high-concentrate diet for 2 weeks 
exhibited a significantly lower ruminal pH, which decreased to 
<5.8 and persisted for more than 3 hours per day after feeding. 
According to the definition of experimental SARA, HC goats were 
suffered SARA disease. However, although ruminal pH value of 
HC goats remained significantly lower than LC goats from week 4 
[4:W) to week 8 {8W), the ruminal pH in HC goats was markedly 
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Figure 3. Hierarchical cluster analysis of the differentially 
expressed transcripts (n = 155). 155 transcripts were differentially 
expressed {adjusted p<0.01). Red and green colors reflect the high and 
low intensities, respectively. 
doi:10.1371/journal.pone.0082850.g003 
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increased to more than 6.0. The increase in ruminal pH of HC 
goats indicated that the rumen was adapting to the high- 
concentrate diet, which was consistent with the study conducted 
in dairy cows by Steele et al [33]. 

Inflammatory response and HPA axis activation 

It is widely accepted that during experimentally-induced SARA, 
the free LPS contents were increased in rumen fluids, yet LPS 
concentration in blood was vague after feeding HC diet that has 
been enriched with starch [15]. A decline in ruminal pH during 
SARA causes death and cell lysis of Gram-negative bacteria, the 
predominant bacterial group in the rumen, which results in an 
increase in free ruminal LPS concentration [12] [15,17]. 
Endotoxin produced in the digestive tract can be translocated 
into the bloodstream, which increases the concentration of blood 
LPS [5] [7]. It's generally accepted that LPS translocation into the 
bloodstream occurs as a consequence of grain-induced SARA, and 
that the intestines are likely the main site of LPS translocation [15] 
[16]. In grain-based SARA, more starch may flow into the hind 
gut especially into the large intestine where LPS may be produced 
in a significant amount as in the rumen via microorganism 
fermentation [9]. After LPS translocates into the peripheral 
circulatory system, the inflammatory and acute phase response will 
be activated [12]. Our results showed that, LPS concentrations 
were increased from 38,854 to 48,064 EU/mL before feeding in 
the morning and increased from 63,359 to 70,418 EU/mL after 
8h feeding compared to LC group. Plasma LPS concentration of 
HC goats was increased from 0.24 to 0.32 EU/mL before feeding 
in the morning and from 0.09 to 0.14 EU/mL after 4 h feeding 
compared with LC group. However, due to the limited animal 
sizes and the individual variations the difference of LPS 
concentrations in ruminal fluids and plasma between HC and 
LC goats did not reach the statistical significance, it can't be 
discarded whether this moderate increase of LPS in body fluids 
involved in activation of inflammatory response. The activation of 
TLR4 and IL-1(3 mRNA expression in liver of HC goats also 
supported the above speculation. Compared to LC group, HC 
goats showed a significant up-regulation in mRNA expression of 
APPs, including hepatic SAA, CRP, and HP mRNA, and 
paralleled a marked increase of SAA and HP in the plasma. 
These results were similar to results reported in dairy cattle [7] 
[13] [15]. 

The previous studies showed that the stress axis is activated 
when SARA is induced [34], Cytokines and APPs are strong 
inducers that signal the HPA axis to release glucocorticoids from 
the adrenal gland [11] [35] [36]. In the present study, we found 
that goats fed a high-concentrate diet for 9 weeks had a 
significandy higher concentration of peripheral blood Cortisol, 
which indicated that the HPA axis was activated during the long- 
term SARA challenge. This is consistent with limited, preliminary 
evidence that high concentrate diet (70% dry matter) adjustment 
to one-year-old Holstein heifers for 30 days increased blood total 
glucocorticoids from 7.9 to 14.8 ng/mL [37]. It's suggested that 
the increased concentration of glucocorticoids may be reflecting a 
metabolic adaptation to high concentrate feed [37]. In this study, 
hepatic GR protein expression was significandy down-regulated in 
liver of HC goats indicating a negative feedback from the higher 
level of plasma Cortisol. As a primary and vital organ that is 
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Figure 4. Co-expression network for Ic-core genes in liver of HC goats. Co-expression networks were for immune and inflammatory 
responses (A), lipid metabolism (B), protein (or amino acid) metabolism (C), and carbohydrate metabolism (D) in the liver of HC goats. Circles in the 
figure represent genes (red circle means up-regulation; blue circle means down-regulation), lines represent regulation relationships between genes 
(solid lines mean a positive correlation of the expression between genes; dotted lines mean a negative correlation), rectangles represent genes for 
immunity and stress, triangles represent carbohydrate metabolism, hexagons represent lipid metabolism, and diamonds represent protein 
metabolism. 

doi:1 0.1 371 /journal.pone.0082850.g004 



involved in the immune response and in nutrient metabolism, the 
liver controls the input of precursors of milk components to the 
mammary gland [24] [38]. These changes of HPA axis may be 
aimed on redistribution of energy utilization in specific organs, 
inhibiting or stimulating energy mobilization according to priority 
[26]. 

The profile of mRNA expression in liver 

In this study, microarray technology was used to investigate the 
hepatic gene expression profile and high-throughput data analysis 
was used to compare the differential expression of functional genes 
between HC and TC group. The Co-expression network of LC 
group shows in Figures SI, S2, S3, and S4. The activation of 
inflammatory and stress response in HC group was also confirmed 
from the GeneRelNet (Co-expression network) dataset. The results 
showed that the significant up-regulation (FC >2, p<0.05) of key- 
core gene regulated genes included GDI (]MM_00 11 23002) and 
CD4 (NM_00 11 29902) molecules, IL-2 pC55641.1) and cortico- 
tropin releasing hormone binding protein (CRHBP, 
NM_001 009339) in liver of HC goats. These genes are important 
for immune response, cell adhesion, and antigen presentation. It's 



well known that CD4 molecule plays a key role for IL-2 produce 
[39] and CDl protein is important for the innate and adaptive 
immune response [40]. Moreover, the significant up-regulation of 
3-hydroxybutyrate dehydrogenase 1 (BDHl, NM_001034600.1) 
and Dipeptidase 1 (DPEPl, NM_001034472.1), and the down- 
regulation of glutathione peroxidase 2 (GPX2, NM_001 163139.1) 
and aldehyde dehydrogenase 18A1 (ALDH18A1, 
NM_001046181.1) indicated the decrease of anti-oxidative 
capacity in liver of goat feeding high concentrate diet [41] [42] 
[43]. In a good agreement, feeding high concentrate diets enrich 
starch to dairy cow induced SARA disease and led to systemic 
inflammatory response due to the translocation of LPS from 
gastrointestinal tract into the peripheral circulation [10] [41] [42]. 

In addition to being an energy source, protein plays a functional 
role in many organs including the liver. The changes in protein 
metabolism are critical to maintaining cell size and directly 
impacts cell survive. Gene ontology analysis also identified 26 
differential expression genes related to protein metabolism. 
Majority of these genes are involved in the biological processes 
of proteolysis, collagen catabolism, and purine nucleotide metab- 
olism. Among these genes, two ATPase (NM_001083653.1, 
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Figure 5. Hepatic gene expression profile detected by real-time 
PCR. (A) mRNA expression of the key genes for inflammatory response 
and acute phase responses. (B) mRNA expression of the key genes for 
glucometabolism. (C) mRNA expression of the key genes for lipid 
metabolism. Beta-actin was used as the reference gene for gene 
expression. Values are mean ± SEM. *p<0.1, *p<0.05 vs. LC. 
doi:1 0.1 371 /journal.pone.0082850.g005 

XM_00 1788345.1), branched chain amino-acid transaminase 1 
(BCATl, NM_00 108 3644.1) and fascin homolog 3 (FSCN3, 
NM_00 1075543.2) were significantly up-regulated in HC goats, 
while RNA polymerase III (NM_001083750.1), and three kinases 
PCM_866451.4, NM_001 109962.1, NM_001078008.1) were 
markedly down-regulated. These changes likely indicated the 
priority of protein catabolism in liver of HC goats compared to 
LC. It's reported that the elevation of plasma Cortisol levels inhibits 
protein synthesis while increasing protein breakdown [43] . In this 
study, the relationship between increased plasma Cortisol level and 
protein breakdown stUl warrants further investigation. 

The priority of catabolism was also observed in lipids 
metabolism. Our results showed that the transcription of genes 
encoding three key enzymes acting in fatty acids P-oxidation [44] , 
ACAA2, ACSLl and HADH were up-regulated in liver of HC 



Table 5. Nutrient composition and forge to concentrate ratio 
(F:C) of the TMR. 




Diet 


Items 


LC 


HC 


LGymus chincnsis 


40 


26 7 


Alfalfa silage 


20 


13.3 


Corn 


22.99 


23.24 


Bran 


0 


20.77 


Soybean meal 


15 


13.66 


Limestone 


0.65 


1.43 


CaHP04 


0.46 


0 


Premix^ 


0.5 


0.5 


NaCI 


0.4 


0.4 


Total 


100 


100 


Nutrient levels'' 


ME (MJ/kg) 


5.63 


5.83 


DCP % 


9.9 


10 


NDF % 


36.64 


34.55 


ADF % 


24.74 


20.35 


EE% 


2.87 


3.21 


NFC % 


31.76 


35 


Ca % 


0.8 


0.9 


P % 


0.33 


0.38 


H20% 


88.4 


88.1 





"Amount of premix added, VA: 6000 lU/kg, VD: 2500 lU/kg, VE: 80 mg/kg, Cu: 
6.25 mg/kg, Fe: 62.5 mg/kg, Zn: 62.5 mg/kg, Mn: 50 mg/kg, I: 0.25 mg/kg, Se: 
0.125 mg/kg, Co: 0.125 mg/kg. Mo: 0.125 mg/kg. 
''Nutrient levels were calculated values. 
LC, low concentrate; HC, high concentrate. 
doi:l 0.1 371 /journa!.pone.0082850.t005 



goats, which was consistent with the increase of NEFA level in 
blood. As a metabolic hormone, leptin plays vital roles on food 
intake, energy expenditure, and body weight mediated by leptin 
receptors situated in the hypothalamus [45]. Leptin has also been 
shown to have effects on whole body glucose and lipid metabohsm 
[46] . In this study, the significant increase of leptin in HC goats 
paralleled the increase of NEFA, which indicates the stimulatory 
effect of leptin on lipids catabolism. With respect to glucose 
metabolism, key-core gene glycogen synthase (GYS2, 
NM_001 192905.1) mRNA transcription was down-regulated in 
HC goats. In addition, Q;PCR results showed that mRNA 
expression of two important gluconeogenic genes hexokinase (HK) 
and pyruvate carboxylase (PC) genes was significantly down- 
regulated in liver of HC group. These data suggests the decrease of 
glucogen synthesis and gluconeogenesis in liver of HC group. 
However, this is contradictory to the increase of blood glucose 
concentrations, which was significantiy elevated in HC group. 
This increase of blood glucose might be related to the production 
of lactose in these lactating animals [47]. Previous study showed 
that fatty acids have shown to produce an insulin-resistant state via 
cascade signaling pathway involving protein kinase C (PKC) and 
AKT phospholation [48]. Also, high levels of free fatty acids in the 
plasma have been shown to induce insulin resistance [49]. 
However, whether the increase of blood glucose in HC group is 
due to insulin-resistance is still unclear. 
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Table 6. PCR primers for immunity and stress, lipid metabolism and carbohydrate metabolism. 





Target genes 


Genbank accession 


PCR products 
(bp) 


Primer sequences 


P-actin 


AF_481159 


260 


F:5'-CGGGATCCATCCTGCGTCTGGACCTG-3' 


R:5 '-GG AATTCGGAAGGAAGGCTGG AAGAG-3 ' 


ACAA2 


XM_004020663.1 


165 


F:5 '-TGTCTGCTGGCAAAGTCTCACC-3 ' 


R:5'-AACCAGAGCCACAGAGCCTGTT-3' 


HADH 


XM_004009637.1 


197 


F:5 '-GAGTTTGTGGCGAAG ACCCTGA-3 ' 


R:5 '-GGCTTGTG ATCTGCAAAGAGG AAG-3 ' 


CPTl 


Y1 8830.1 


158 


F:5'-GCTTTGATCGACACTTGTTTGCTC-3' 


R:5 ■ -GC A A AGC AGCCG ATGTTC ACT-3 ' 


CPT2 


BC 105423.1 


157 


F:5'-GCTTTGATCGACACTTGTTTGCTC-3' 


R:5 ■-GCAAAGCAGCCGATGTTCACT-3 ' 


ACSLl 


NM_001 076085.1 


191 


F:5 '-GCAACCCCAAAGGAGCAATG-3 ' 


R:5'-AGCAGCCTGATATCTCCTTGG-3' 


PPARa 


HM600810.1 


243 


F:5 '-TAAAGCCAACCAAG ATAACCC-3 ' 


R:5 '-TCACCAAACAGCCGA AGA-3 ' 


ATGL 


GQ918145.1 


180 


F:5'-GGAGCTTATCCAGGCCAAT-3' 


R:5'-TGCGGGCAGATGTCACTC-3' 


ACSS1 


BC055008.1 


242 


F:5 '-TCCTTGGCTGGG AGGATCAA-3 ' 


R:5 '-TGTTGTCTGTCCTGTG AGCCA-3 ' 


LPL 


JQ670882.1 


235 


F:5'-TTCAGAGGCTATTACTGGAAATCC-3' 


R:5 ■ -ATGTC A ATC AC AGC ATTC ATTCTACT-3 ' 


SCD 


AF325499.1 


178 


F:5'-TGCTGACAACTTATCTGGATGC-3' 


R:5'-AAGGAATCCTGCAAACAGCTA-3' 


DGATl 


DQ380249.1 


240 


F:5 ' -TGCCTC AG AC ACTTCTAC A AGCC-3 ' 


R:5 '-GCCCG ATG ATG AGTG ACAGC-3 ' 


DGAT2 


AJ51 9787.1 


234 


F:5 '-CACTGGCTCCAGCATCCTCTC-3 ' 


R:5 '-TTCTTGGGTGTGTTCCAGTCA-3 ' 


FADSl 


XM_004019593.1 


187 


F:5 '-CCTTGCTGCCTCTATACTTCCA-3 ' 


R:5 '-ACAAACC AGTTGCTTTCCAGGA-3 ' 


FADS2 


AY.850395.1 


257 


F:5'-AGAGCATCGCCTGGTTCACTA-3' 


R:5 '-CCTTGTGG AAG ACATTGGGTTT-3 ' 


ACACA 


NM_001 009256 


230 


F:5'-CATGGAAATGTACGCGGACC-3' 


R:5 '-GGTGGTAGATGGG AAGGAGGA-3 ' 


FAS 


XM_00401 3447.1 


112 


F:5 '-TGCTCATTCACTCGGGTTCT-3 ' 


R:5 ■ -AGGTATGCCCGCTTTTCG-3 ' 


PPARy 


JN854246.1 


238 


F:5'-CATTTCTGCTCCGCACTAC-3' 


R:5 '-TGGAACCCTGACGCTTT-3 ' 


PC 


NM_1 77946 


187 


F:5'-TCGCACCATGTATGTCATCCC-3' 


R:5'-AGGCI II II 1 AAAGGCAGAGGG-3' 


FBPl 


NM_00 1034447 


106 


F:5'-GCGGTCAAAGCCATCTCCAC-3' 


R:5 '-CATCCAGCTTCTTCACTTG ATCTCC-3 ' 


G6Pase 


EF062861 


158 


F:5'-AATGTCATGTTGTGGTTGGGA"n"CT-3' 


R:5 '-GCATTGTAG ATGCTCTGG ATGTGG-3 ' 


HK1 


AM492192 


139 


F:5 '-GCGGCTCTCTGATAAAACTCTGTTA-3 ' 


R:5 '-TGAGCCATCGGG AATAG ACCTTAC-3 ' 


PYGL 


AY827551 


150 


F:5 ' -GCCTTCCC AG ATC AGGTTGC-3 ' 


R:5'-GTGGTTGGTGTAGGCGAAGGTC-3' 


LEPR 


NM_001 009763 


110 


F:5'-GGAAGGAGTAGGGAAACCGAAGA-3' 


R:5 '-TTGAGGAGG AG ATTATTATTGGCAC-3 ' 


GYS2 


NM_001 192905 


131 


F:5'-TTGGGCGGTATCTTTGTG-3' 
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Table 6. Cont. 



PCR products 

Target genes Genbank accession (bp) Primer sequences 

R:5'-AGCACAATGGACGGAAGC-3' 
GSKSa NM_001 102192.1 250 F:5'-TGGCTTACACAGACATCAAA-3' 

R:5'-TCGGGCACATATTCCAGCAC-3' 
GSK3|3 NM_001 101310.1 249 F:5'-AGACAAAGATGGCAGCAAGGTGAC-3' 

R:5 '-ACGCAATCGGACTATGITAC-S ' 

GLUT2 NM_001 103222.1 268 F:5'-CATCCATC"n"CCTC"nTGTCTG-3' 

F:5 '-GATTTTCCTTTGGTTTCTGG-3 ' 
IR XM_590552.4 275 F:5'-GACGCAGGCCGGAGATGACCA-3' 

R:5 '-GCTCCTGCCCGAAG ACCG ACTC-3 ' 
SAA AF540564.1 121 F:5'-CATCCTGCGTCTGGACCTGG-3' 

R:5 '-TTCCTTG ATGTCACGGACG ATTT-3 ' 
TLR4 JQ342090.1 195 F: 5'-G"nTCCACAAGAGCCGTAA-3' 

R: 5'- TGTTCAGAAGGCGATAGAGT-3' 

MyD88 JQ308783.1 98 F:5'-ACAAGCCAATGAAGAAAGAG-3' 

R:5 '-G AGGCGAGTCCAG AACC-3 ' 
CD14 NM_001 077209.1 239 F:5'-CCGTTCAGTGTATGGTTGCC-3' 

R:5 ■-TGCTTCGGGTCGGTGTT-3 ' 
HP XM_004015111.1 162 F:5'-TAATGCCCATCTGCCTAC-3' 

R:5 ■ -CGCCCTC ATAGTGTTTC A-3 ' 
IL-ip D63351.1 172 F: 5'-GAAGAGCTGCACCCAACA-3' 

R: 5'- CAGGTCATCATCACGGAAG-3' 
TNF-a AF276985.1 173 F: 5'-CAAGTAACAAGCCGGTAGCCC-3' 

R:5'-CCTGAAGAGGACCTGCGAGTAG-3' 
LBP XM_004014566.1 138 F: 5'-TGGAGCCAGGAAAGATAC-3' 

F: 5'-AGCCTTCTGCCAACTTAT-3' 
CRP NM_00 1144097.1 134 F: 5'-CTGGCTTGGGAGATTG-3' 

F: 5'-AGTGAGGGTAAGGGATT-3' 

ACAA2, acetyl-coenzyme A acyltransferase 2; ACACA, acetyl-coenzyme A carboxylase alpha; ATGL, adipose triglyceride lipase; ACSLl, acyl-CoA synthetase long-chain 
family member 1; ACSSl, acyl-CoA synthetase short-chain family member 1; CD14, CD14 molecule; CPT1, carnitine palmltyl transferasel; CPT2, carnitine palmityl 
transferase 2; CRP, C-reactive protein; DGATl, diacylglycerol acyltransferase 1; DGAT2, diacylglycerol acyltransferase 2; FADSl, fatty acid desaturase 1; FADS2, fatty acid 
desaturase 2; FBP1, fructose-l,6-bisphosphatase 1; PPARa, peroxisome proliferator-activated receptor a; FAS, fatty acid synthase; G6Pase, glucose-6-phosphatase; 
GLUT2, solute carrier family 2 {facilitated glucose transporter), member 2; GYS2, glycogen synthase 2; GYSK3q(, glycogen synthase kinase 3 alpha; GYSK3P, glycogen 
synthase kinase 3 beta; HADH, 3-hydroxyacyl-CoA dehydrogenases; HKl, hexokinase 1; HP, haptoglobin; ILIP; interleukin-lbeta; IR, insulin receptor; LBP, 
lipopolysaccharide binding protein; LEPR, leptin receptor; LPL, lipoprotein lipase; MyD88, myd88; PC, pyruvate carboxylase; PPARy, peroxisome proliferators-activated 
receptor y; PYGL, liver glycogen phosphorylase; SAA, serum amyloid A protein; SCD, stearoyl-CoA desaturase; TLR4, toll-like receptor 4; TNFa, tumor necrosis factor 
alpha. 

doi:1 0.1 371 /journal.pone.0082850.t006 



Taken as a whole, our data suggest that feeding HC diet to 
lactating goats for long period caused the activation of the 
inflammatory and stress response, and altered substrates metab- 
olism in liver. Based on the data presented in this study, limited 
conclusions can be drawn from just measuring mRNA abundance 
because the post-transcription, protein translation, and the 
amount of active protein are dependent on many factors [50]. 
To our knowledge, it's the first time to report the changes of 
metabolism in Uver after feeding HC diet to lactating ruminants, 
and the determination of wide-genome expression in liver will 
contribute to find the specific biochemical pathways. 

Conclusion 

In summary, feeding a high-concentrate diet to lactating goats 
for the long-term will result in the depression of milk fat and milk 
protein. Activation of the stress response, the acute phase response 



stimulation and the inflammatory response may contribute to 
nutrient partitioning and re-distribution of energy in the liver, and 
ultimately lead to a decline in milk quality. 

Materials and Methods 

Ethics Statement 

The experiment was conducted following the guidelines of 
Animal Ethics Committee at Nanjing Agricultural University, 
China. The study was approved by Animal Ethics of Nanjing 
Agriculture University. The sampling procedures complied with 
the the "Guidelines on Ethical Treatment of Experimental 
Animals" (2006) No. 398 set by the Ministry of Science and 
Technology, China and "the Regulation regarding the Manage- 
ment and Treatment of Experimental Animals" (2008) No. 45 set 
by the Jiangsu Provincial People's Government. 
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Animals 

Eight health multiparous mid-locating goats with average initial 
BW 47 ±8 kg (mean ± SD) were housed in individual stalls in a 
standard animal feeding house at Nanjing Agricultural University 
(Nanjing, China). Animals were cannulated in rumen and were 
randomly allocated to two groups, one receiving diets with low 
concentrate (40% of dry matter) as the LC group (n = 4), and 
another receiving high concentrate diet (60% of dry matter) as the 
HG group (n = 4) for 9 wcc'ks. Ingredient and nutrient composition 
of experimental diets were presented in Table 5. Goats were free 
access to fresh water throughout the experimental time. 

Milk composition analysis 

Fifty milliliter of milk were sampled every two weeks into the 
vials with potassium dichromate at 8:00 in the morning, and stored 
at the 4°C until milk fat, protein and lactose analysis conducted by 
the commercied company of Nanjing Weigang dairy industry Co., 
Ltd. The data were analyzed by Univariate using the General 
Linear. Models of SPSS 11.0 for Windows (StatSoft Inc, Tulsa, 
OK, USA). 

Rumen fluid collecting and measurements 

Fifteen minutes prior to feed delivery and 1, 2, 3 and 5 h (or 1, 
3, 5 h) after ftx-d d(4ivery on 2 consecutive days during 2, 4, 6 and 
8 week of experimental period, 15 mL rumen fluids was collected 
with nylon bag and divided into 5 mL for pH measurement 
immediately with pH-meter. Another 1 0 mL were transferred into 
tubes and kept on ice and then centrifuged at 1,900 xg at 4°C for 
30 min. After that, the supernatant was collected and stored at 
— 20°C until for LPS analysis. The concentration of LPS in rumen 
fluid was determined by a chromogenic endpoint assay (CE64406, 
Chinese Horseshoe Crab Reagent Manufactory Co., Ltd., Xia- 
men, China) with a minimum detection limit of 0. 1 EU/ mL. The 
procedures were performed according to the manufacture's 
instruction. The data of pH and LPS concentrations were 
analyzed by Univariate using the General Linear Models of 
SPASS 11.0 for windows (StatSoft Inc, Tulsa, OK, USA). 

Measurement of plasma biochemical parameters 

At the end of experiment, plasma was sampled fifteen minutes 
prior to feed delivery using EDTA-containing vacuum tubes from 
jugular vein. Plasma leptin concentration was measured with 
radioimmunoassay (RIA) using commercial kits purchased from 
Shanghai Institute of Biological Products. The ranges of assay- 
sensitivity were between 0.5 and 24 ng/mL. The inter- and intra- 
assay coefficients of variation were 5% and 10%, respecti\(-ly. Tlu- 
concentration of NEFA (A042) and glucose was detected by 
commercial kits purchased from Nanjing Jiancheng Bioengineer- 
ing institution, and the procedures were performed according to 
the manufacture's instruction. The concentrations of leptin and 
glucose were carried out in Nanjing General Hospital of Nanjing 
Military. Inflammatory cytokines such as IL-6 (RGB-60023G, 
Beijing Rigorbio Science Development Co., Ltd.), TNFa (RGB- 
60080G) and acute phase protein such as LBP (RGB-60266G), HP 
(H086-09, abeam), SAA (ab 100635, abeam) were detected by 
ELISA kits. The concentration of LPS in plasma was determined 
by a chromogenic endpoint assay with diazo coupling reagent 
(CE80545, Chinese Horseshoe Crab Reagent Manufactory Co., 
Ltd., Xiamen, China) with a minimum detection Kmit of 
0.01 EU/mL. The procedures were performed according to the 
manufacture's instruction. The data were analyzed by Indepen- 
dent-Samples T test using the Compare Means of SPASS 1 1.0 for 
Windows (StaSoft Inc, Tulsa, OK, USA). 



RNA extraction and real-time quantitative PCR 

The liver tissues were sampled by a punch biopsy with a local 
anesthesia. Total RNA was extracted from each dairy goat using 
TRIZOL reagent (Takara, Dalian, China) according to the 
manufacturer's protocols. Total RNA concentration was then 
quantified by measuring the absorbance at 260 nm in a 
spectrophotometer (Eppendorf Biotechnology, Hamburg, Ger- 
many). Aliquots of RNA samples were subjected to electrophoresis 
with 1.4% agarose-formaldehyde gels stained with ethidium 
bromide to verify their integrity. Reverse transcription (RT) was 
performed using the total RNA (2 |xg) in a final volume of 25 |tL 
containing 1 x RT-buffer, 100 U reverse transcriptase, 8 U RNase 
inhibitor (Promega, USA), 5.3 (imol/L random hexamer primers 
and 0.8 mmol/L dNTP (TaKaRa, Dalian, China). After incuba- 
tion at 37°C for 1 h, the reaction was terminated by heating at 
95°C for 5 min and quickly cooling on ice. Real-time PCR was 
performed in Mx3000P (Stratagene, La JoUa, CA, USA). Mock 
RT and No Template Controls (NTC) were set to monitor the 
possible contamination of genomic DNA both at RT and PCR. 
Melting curves were performed to insure a single specific PCR 
product for each gene. Two microliter of 16-fold dilution of RT 
product was used for PCR in a final volume of 25 |iL containing 
12.5 [iL SYBR Green Real-time PCR Master Mix (TOYOBO 
Ltd., Shanghai, China) and 0.6-0.8 [lM of each forward and 
reverse primers for target genes (as shown in Table 6). Goat P- 
actin mRNA was used as a reference gene for normalization 
purposes. The following PCR protocols were initial denaturation 
(1 min at 95°C), then a three-step amplification program (20 s at 
95°C, 20-30 s at 60-64°C, 30 s at 72°C) was repeated 45 times. 
The method of 2 ^^'^^ was used to analyze qPCR data with 
Independent-Samples T test using Compare means of SPSS 11.0 
for Windows to statistical analysis. AU samples were included in 
the same run of RT-PCR and repeated for at least 3 times. The 
primers were shown in Table 6. The data was analyzed by 
Independent-Samples T test using the Compare Means of SPASS 
11.0 for Windows (StaSoft Inc, Tulsa, OK, USA). 

Western blotting 

Total protein was extracted from frozen liver tissue, and the 
protein concentration was determined by the BCA assay (Pierce, 
Rockford, IL, USA). Fifty micrograms of protein extract from each 
sample were subjected to electrophoresis on a 10% SDS-PAGE 
gel, and the separated proteins were transferred onto the 
nitrocellulose membranes (Bio Trace, Pall Co., USA). Western 
blot analysis for GR (sc-1004, Santa Cruz Biotechnology, 1:500) is 
performed with primary antibodies and corresponding HRP- 
conjugated secondary antibodies. P-actin (KC-5A08, KangChen 
Bio-tech, China, 1:5000) was used as a reference protein for 
normalization purposes in the Western blot analysis. Finally, the 
blot was washed and detected by enhanced chemiluminescence 
(ECL) using the LumiGlo substrate (Super Signal West Pico Trial 
Kit, Pierce, USA). ECL signals were recorded by an imaging 
system (Bio-Rad, USA) and analyzed with Quantity One software 
(Bio-Rad, USA). Values of GR protein was presented as fold 
change relative to the average value of LC group. The data was 
analyzed by Independent-Samples T test using the Compare 
Means of SPASS 11.0 for Windows (StaSoft Inc, Tulsa, OK, 
USA). 

mRNA Microarray Experiment 

Transcriptome microarrays were performed for seven dairy 
goats (3 goats from the LC group and 4 goats from the HC group) 
for hepatic gene expression profiling. Because complete genome 
sequence data was not available, gene sequences from three 
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animal species were employed to design the gene chip (goat, sheep, 
and Bos Taurus). Goat expressed sequence tag (EST) sequence 
(n= 13898), goat gene sequence (n = 720), sheep gene sequence 
(n = 2718), and Bos Taurus gene sequence (n> 35000) transcripts 
were investigated using the microarray. There were 61012 probe 
IDs and 2-fold difference probe IDs (n = 2140). Reference 
sequences were cited from httpV/www.ncbi.rdm.nih.gov/gene/ 
and EST database was from http://www.ncbi.nlm.nih.gov/ 
nucest/. The FDR (False Discovery Rate) were calculated. 
Differentially expressed genes (DEGs) were selected with FDR 
<5% and FDR <10%. All data were MIAME compliant and 
have been deposited in GEO (accession number GPL17623 and 
GPL17624). 

GO analysis 

GO analysis was applied to analyze the main function of the 
differential expression genes according to the Gene Ontology 
which is the key functional classification of NCBI, which can 
organize genes into hierarchical categories and uncover the gene 
regulatory network on the basis of biological process and 
molecular function [51]. Specifically, two-side Fisher's exact test 
and test were used to classify the GO category, and the false 
discovery rate (FDR) [52] was calculated to correct the P-value,the 
smaller the FDR, the small the error in judging the p-value. The 

FDR was defined as FDR =1 — — , where Nk refers to the 

number of Fisher's test P-values less than test P-values. We 
computed P-values for the GOs of aU the differential genes. 
Enrichment provides a measure of the significance of the function: 
as the enrichment increases, the corresponding function is more 
specific, which helps us to find those GOs with more concrete 
function description in the experiment. Within the significant 
category, the enrichment Re was given by: ^e = (nf/n)/(Nf/N) 
where is the number of flagged genes within the particular 

category, "«" is the total number of genes within the same 
category, "Nf" is the number of flagged genes in the entire 
microarray, and "N" is the total number of genes in the 
microarray [53]. 



References 

1. Yang VVZ, Brauchcmin KA (2006) Effects of physically effective fiber on 
chewing acli\ity and mminal pH of dairy cows fed diets based on barley silage. 
J Dairy Sri 89: 217-228. 

2. Kcuncn JE, Plaizicr JC, Kyriazakis L, Dufiicld TF, Widowski TM, ct al. (2002) 
Effects of a subacute ruminal acidosis modei on the diet selection of dairy cows. 
J Dairy Sci 85: 3304-3.S1.S. 

3. Kleen JL, Hooijcr GA, Rehage J, Noordhuizen JP (2003) Subacute ruminai 
acidosis (SARA): a review. J Vet Med A Physiol Pathol Chn Med 50: 406-414. 

4. Piaizier JC, Krause DO, Gozho ON, McBride BW (2008) Subacute ruminal 
acidosis in dairy cows: the physiologicai causes, incidence and consequences. 
Vet J 176: 21-31. 

5. Gozho GN, Piaizier JC, Krause DO, Kennedy AD, Wittenberg KM (2005) 
Subacute ruminal acidosis induces ruminal lipopolysaccharide endotoxin release 
and triggers an inflammatory response. J Dairy Sci 88: 1399-1403. 

6. Beauchemin Kj\, Yang \VZ, Rode LM (2003) Effects of particle size of alfalfa- 
based dairy cow diets on chewing activity, ruminal fermentation, and milk 
production. J Dairy' Sci 86: 630—643. 

7. Fairfield AM, Plaizicr JC, Dufiicld TF, Lindingcr Ml, Bagg R, et al. (2007) 
Effects of prepartum administration of a monensin controlled release capsule on 
rumen pH, feed intake, and milk production of transition dairy cows. J Dairy Sci 
90: 937-945. 

8. Khafipour E, Krause DO, Piaizier JC (2009) Alfalfa peUet-induced subacute 
ruminal acidosis in dairy cows increases bacterial endotoxin in the rumen 
without causing inflammation. J Dairy Sci 92: 1712-1724. 



GeneRelNet (Co-expression networl<) 

Gene co-expression Networks was estabfished to identify the 
interactions among genes. Gene co-expression Networks were 
buUt according to the normalized signal inttnisity of specific 
expression genes [54]. For each pair of genes, we calculate the 
Pearson correlation and choose the significant correlation pairs 
with which to construct the network [55]. In a network analysis, 
degree centrality is the simplest and most important measures of a 
gene centrality within a network that determining the relative 
importance. Degree centrality is defined as the link numbers one 
node has to the other [56]. Moreover, to study a variety of 
properties of networks, k-cores were introduced in graph theory as 
a method of simplifying graph topology analysis. A k-core of a 
network is a subnetwork in which all nodes are connected to at 
least k other genes in the subnetwork. The Purpose of Network 
Structure Analysis is to locate core regulatory factors (genes). In 
one network, core regulatory factors connect most adjacent genes 
and have the biggest degrees. While considering different 
networks. Core regulatory factors were determined by the degree 
differences between two class samples [57]. 

Supporting Information 

Figure SI Co-expression network for immune and 
inflammatory responses in the liver of LC goats. 

(TIF) 

Figure S2 Co-expression network for lipid metabolism 
in the liver of LC goats. 

(TIF) 

Figure S3 Co-expression network for protein (or amino 
acid) metabolism in the liver of LC goats. 

(TIF) 

Figure S4 Co-expression network for carbohydrate 
metabolism in the liver of LC goats. 

(TIF) 

Author Contributions 

Conceived and designed the experiments: YN YZ SZ XS. Performed the 
experiments: HD SW YJ. Contributed reagents/materials/analysis tools: 
RZ. Wrote the paper: HD YN. 



9. Li S. Khafipoiir I'j. Kraiisc DC). Krorkrr A, Rodrig'iicz-Lrcomptc JC, ct al. 
(^(11 2': ILIliTl.s dl' Mil)at'ulc ruminal acidosis rliallrngcs on fermentation and 
endotoxins in tlie rumen and hindgut of dairy cows. J Dairy Sci 95: 294-303. 

10. Khafipour E, Krause DO, Piaizier JC (2009) A grain-based subacute ruminal 
acidosis challenge causes translocation of lipopolysaccharide and triggers 
inflammation. J Dairy Sci 92: 1060-1070. 

11. NocekJE (1997) Bovine acidosis: implications on laminitis. J Dairy Sci 80: 1005- 
1028. 

12. Dong G, liu S, Wu Y, Lei C, Zhou J, et al. (2011) Diet-induced bacterial 
immunogens in the gastrointestinal tract of dairy cows: impacts on immunity and 
metabolism. Acta Vet Scand 53; -IS. 

13. Cruys E, Toussaint MJ, Niewold TA, Koopmans SJ (2005) Acute phase reaction 
and acute phase proteins. J Zhejiang Univ Sci B 6: 1045-1056. 

14. Bannerman DD, Paape MJ, Hare WR, Hope JC (2004) Characterization of the 
bovine innate immune response to intramammary infection with Klebsiella 
pneumoniae. J Dairy Sci 87: 2420-2432. 

15. Sohn MJ, Hut GM, Byun HS, Kim WG (2008) Cyclo(dehydrohistidyl-l- 
tryptophyl) inhibits nitric oxide production by preventing the dimerization of 
inducible nitric oxide synthase. Biochem Pharmacol 75: 923-930. 

16. Baumaim H, Gauldie J (1994) The acute phase response. Immunol Today 15: 
74-80. 

17. Heinrich PC, Castell JV, Andus T (1990) Interleuldn-6 and the acute phase 
response. BiochemJ 265: 621—636. 

18. Taguchi Y, Yamamoto M, Yamate T, Lin SC, Mocharla H, et al. (1998) 
Interleukin-6-type cytokines stimulate mesenchymal progenitor differentiation 
toward the osteoblastic lineage. Proc Assoc Am Physicians 1 10: 559-574. 



PLOS ONE I www.plosone.org 



11 



December 2013 | Volume 8 | Issue 12 | e82850 



Long-Term Effects of SARA on Milk Quality 



19. Murata H, Shimada N, Yoshioka M (2004) Current research on acute phase 
proteins in veterinary diagnosis: an overview. Vet J 168: 28—40. 

20. Heinrich PC, Behrmann I, MuUer-Newen G, Schaper F, Graeve L (1998) 
Interleukin-6-type cytokine signalling through the gpl30/Jak/STAT pathway. 
Biochem J 334 (Pt 2): 297-314. 

21. Zunszain PA, Anacker C, Cattaneo A, Carvalho LA, Pariante CM (2011) 
Glucocorticoids, cytokines and brain abnormalities in depression. Prog 
Neuropsychopharmacol Biol Psychiatry 35: 722—729. 

22. Dentin R, Girard J, Postic C (2005) Carbohydrate responsive element binding 
protein (ChREBP) and sterol regulatory element binding protein- Ic (SREBP-lc): 
two key regulators of glucose metabolism and lipid synthesis in liver. Biochimie 
87: 81-86.^ 

23. Johnson RW (1997) Inhibition of growth by pro-inflammatory cytokines: an 
integrated view. J Anim Sci 75: 1244-1255. 

24. Nukina H, Sudo N, Aiba Y, Oyama N, Koga Y, et al. (2001) Restraint stress 
elevates the plasma interleuldn-6 levels in germ-free mice. J Neuroimmunol 1 15: 
46-52. 

25. Alscmgcest SP, Lambooy IE, Wierenga HK, Dielcman SJ, Meerkerk B, et al. 
(1995) Influence of physical stress on the plasma concentration of serum 
amyloid-A (SAA) and haptoglobin (Hp) in calves. Vet Q_ 17: 9-12. 

26. C;rum DE, Drackley JK, Hansen LR, Cremin JD Jr. (1996) Production, 
digestion, and hepatic lipid metabolism of dairy cows fed increased energy from 
fat or concentrate. J Dairy Sci 79: 1836-1849. 

27. Steele MA, Alzahal O, Walpole ME, McBride BW (2012) Short communication: 
grain-induced subacute ruminal acidosis is associated with the differential 
expression of insulin-like growth factor-binding proteins in rumen papiUae of 
lactating dairy cattle. J Dairy Sci 95: 6072-6076. 

28. Garrett EE, Pereira MN, Nordlund KV, Armentano EE, (Joodger WJ, ct al. 
(1999) Diagnostic methods for the delertion ol'siibaeiile ruminal acidosis in dairy 
cows. J Dairy Sci 82: 1 1 70-1 1 78. 

29. Oba M, Allen MS (2000) Eflcels of brown midrib 3 mutation in corn silage on 
productivity of dairy cows fed u\o concentrations of dietary neutral detergent 
fiber: 3. Digestibility and microbial cflicicncy. J Dairy Sci 83: 1350-1358. 

30. AlZahal O, Rustomo B, Odongo NE, Duffield TF, McBride BW (2007) 
Technical note: A system for continuous recording of ruminal pH in catde. 
J Anim Sci 85: 213-217. 

31. DufTield T, Plaizier JC, Fairfield A, Bagg R, Vessie G, et al. (2004) Comparison 
of techniques for measurement of rumen pH in lactating dairy cows. J Dairy Sci 
87: 59-66. 

32. Yang \VZ, Beauchemin KA (2006) Effects of physically effective fiber on 
chewing activity and ruminal pH of dairy cows fed diets based on barley silage. 
J Dairy Sci 89: 217-228. 

33. Steele MA, Croom J, Kahler M, AlZahal O, Hook SE, et al. (2011) Bovine 
rumen epithelium undergoes rapid structural adaptations during grain-induced 
subacute ruminEil acidosis. Am J Physiol Regul Integr Comp Physiol 300: 
R1515-1523. 

34. Minton JE (1994) Function of the hypothalamic-pituitary-adrenal axis and the 
sympathetic ner^'ous system in models of acute stress in domestic farm animals. 
J Anim Sci 72: 1891-1898. 

35. Kabaroff LC, Rodriguez A, Quinton M, Boermans H, Karrow NA (2006) 
Assessment of the ovine acute phase response and hepatic gene expression in 
response to Escherichia coli endotoxin. Vet Immunol Immunopathol 113: 113— 
124. 

36. Karrow NA (2006) Activation of the hypothalamic-pituitary-adrenal axis and 
autonomic nervous system during inflammation and altered programming of the 
neuroendocrine-immune axis during fetal and neonatal development: lessons 
learned from the model inflammagen, Hpopolysaccharide. Brain Behav Immun 
20: 144-158. 



37. Mills SE, Jenny BE (1979) Eliects of high concentrate feeding and fasting on 
plasma glucocorticoids in dairy heifers. J Anim Sci 48: 961—965. 

38. Sander EE, Sackett SD, Dierssen U, Beraza N, linke RP, et al. (2010) Hepatic 
acute-phase proteins control innate immune responses during infection by 
promoting myeloid-derived suppressor cell function. J Exp Med 207: 1453— 
1464. 

39. Zhu J, Yamane H, Paul WE (2010) Differentiation of effector CD4 T cell 
populations (*). Annu Rev Immunol 28: 445-489. 

40. Hussam MM, Rava P, Walsh M, Rana M, Iqbal J (2012) Multiple functions of 
microsomal triglyceride transfer protein. Nutr Adetab (Lond) 9: 14. 

41. Emmanuel DG, Dunn SM, Ametaj BN (2008) Feeding high proportions of 
barley grain stimulates an inflammatory response in dairy cows. J Dairy Sci 91: 
606-614. 

42. GN G, Plaizier Jc Fan - Krause DO, Krause Do Fau - Kennedy AD. Kennedy 
Ad Fau - Wittenberg KM, Wittenberg KM (2005) Subacute ruminal acidosis 
induces ruminal hpopolysaccharide endotoxin release and triggers an inflam- 
matory response. 

43. Chen (iL, Miller GM (2012) Advances in tryptophan hydroxvlasc-2 gene 
expression regulation: new insights into serotonin-stress interaction and clinical 
implications. AmJ Med Genet B Neuropsychiatr Genet 159B: 152-171. 

44. Mashek DG, Li LO, Coleman RA (2007) Long-chain acyl-CoA synthetases and 
fatty acid channeling. Future Lipidol 2: 465—476. 

45. Stephens TW, Basinski M, Bristow PK, Bue-Valleskey JM, Burgett SG, et al. 
(1995) The role of neuropeptide Y in the antiobesity action of the obese gene 
product. Nature 377: 530-532. 

46. BrysonJM, PhuyalJL, Swan V, Caterson ID (1999) Leptin has acute effects on 
glucose and lipid metabolism in both lean and gold thioglucosc-obcse mice. 
AmJ Physiol 277: E417-422. 

47. Xiao CT, Cant JP (2005) Relationship between glucose transport and 
metabolism in isolated bovine mEunmary epithelial cells. J DEury Sci 88: 2794— 
2805. 

48. Gao Z, Zhang X, Zuberi A, Hwang D, Quon A4J, et al. (2004) Inhibition of 
insulin sensitivity by free fatty acids requires activation of multiple serine kinases 
in 3T3-L1 adipocytes. Mol Endoerniol 18: 2024-2034. 

49. Boden (i, Homko G, Mozzoli M, Zhang M, Kresge K, et al. (2007) Combined 
use of rosiglitazone and fenofibratc in patients with t\'jie 2 diabetes: prevention of 
fluid retention. Diabetes 56: 248-255. 

50. Salisbury J, Hutchison KW, Wigglesworth K, EppigJJ, Graber JH (2009) Probe- 
level analysis of expression microarrays characterizes isoform-specific degrada- 
tion during mouse oocyte maturation. PLoS One 4: e7479. 

51. Ashbumer M, Ball CA, Blake JA, Botstein D, Butler H, et al. (2000) Gene 
ontology: tool for the unification of biology. The Gene Ontology Consortium. 
Nat Genet 25: 25-29. 

52. Dupuy D, Bcrtin N, Hidalgo CA, Venkatesan K, Tu D, et al. (2007) (Genome- 
scale analysis of in vivo spatiotemporal promoter activity in Gacnorhabditis 
elegans. Nat Biotcchnol 25: 663-668. 

53. Schlitt T, Palin K, Rung J, Dietmann S, Lappc M, et al. (2003) From gene 
networks to gene function. Genome Res 13: 2568-2576. 

54. Pujana MA, Han JD, Starita LM, Stevens KN, Tewari M, et al. (2007) Network 
modeling links breast cancer susceptibility and centrosome dysfunction. Nat 
Genet 39: 1338-1349. 

55. Prieto C, Risueno A, Fontanillo C, De las Rivas J (2008) Human gene 
coexpression landscape: confident network derived from tissue transcriptomic 
profiles. PLoS One 3: e3911. 

56. Barabasi AL, Oltvai ZN (2004) Network biology: understanding the cell's 
functional organization. Nat Rev Genet 5: 101-113. 

57. Carlson MR, Zhang B, Fang Z, Mischel PS, Horvath S, et al. (2006) Gene 
connectivity, function, and sequence conservation: predictions from modular 
yeast co-expression networks. BMC Genomics 7: 40. 



PLOS ONE I www.plosone.org 



12 



December 2013 | Volume 8 | Issue 12 | e82850 



